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Abstract: We consider the problem of reconstructing binary images from their horizontal and 
vertical projections. For any reconstruction we define the length of the boundary of the image. In 
this paper we assume that the projections are monotone, and we construct an image satisfying 
these projections that has a relatively small boundary. We also give families of examples for which 
we show that no smaller boundary is possible. 



1 Introduction 

An important problem in discrete tomography is to reconstruct a binary image on a lattice 
from given projections in lattice directions [HI HO] . Each point of a binary image has a value 
equal to zero or one. The line sum of a line through the image is the sum of the values of 
the points on this line. The projection of the image in a certain direction consists of all the 
line sums of the lines through the image in this direction. Any binary image with exactly 
the same projections as the original image we call a reconstruction of the image. 

For any set of more than two directions, the problem of reconstructing a binary image 
from its projections in those directions is NP-complete [7]. For exactly two directions, the 
horizontal and vertical ones, say, it is possible to reconstruct an image in polynomial time. 
Already in 1957, Ryser described an algorithm to do so [12]. He also characterised the set 
of projections that correspond to a unique binary image. 

If there are multiple images corresponding to a given set of line sums, it is interesting 
to reconstruct an image with a special property. In order to find reconstructions that 
look rather like a real object, two special properties in particular are often imposed on the 
reconstructions. The first is connectivity of the points with value one in the picture [2|[5|[T5]. 
The second is hv-convexity: if in each row and each column, the points with value one form 
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one connected block, the image is called hv-convex. The reconstruction of hv-convex images, 
either connected or not necessarily connected, has been studied extensively IH 121 [3l lU [13]. 

Another relevant concept in this context is the boundary of a binary image. The boundary 
can be defined as the set of pairs consisting of two adjacent points, one with value and 
one with value 1. Here we use 4-adjacency: that is, a point is adjacent to its two vertical 
and to its two horizontal neighbours |TT]. The number of such pairs of adjacent points with 
two different values is called the length of the boundary or sometimes the perimeter length 

In this paper we will consider given line sums that may correspond to more than one binary 
image. Since the boundary of real objects is often small compared to the area, it makes 
sense to look for reconstructions of which the length of the boundary is as small as possible. 
In particular, if there exists an hv-convex reconstruction, then the length of the boundary 
of that image is the smallest possible. In that sense, the length of the boundary is a more 
general concept than hv-convexity. 

In [6] we proved a lower bound on the length of the boundary for any reconstruction 
of a picture with given line sums. In this paper we complement this result by giving a 
reconstruction that has a relatively small boundary in the case that both the row and the 
column sums are monotone. 

After introducing some notation in Section |2j we describe the construction of a solution 
to the discrete tomography problem in Section [3j In Section |4] we prove upper bounds on 
the length of the boundary of this constructed solution. We show by examples that these 
bounds are sharp in Section |5| and finally in Section |6] we generalise the results slightly. 



2 Definitions and notation 

Let F be a finite subset of with characteristic function x- (That is, x{k, I) = 1 ii {k, I) G 
F and x{k, I) = otherwise.) For z G Z, we define row i as the set {{k,l) ^T? : k = i}. We 
call i the index of the row. For j G Z, we define column j as the set {{k,l) G Z^ : / = j}. 
We call j the index of the column. Note that we follow matrix notation: we indicate a 
point by first its row index i and then its column index j. Also, we use row numbers 
that increase when going downwards and column numbers that increase when going to the 
right. 

The row sum ri is the number of elements of F in row i, that is = X^jg^xl^jj)- The 
column sum cj of F is the number of elements of F in column j, that is Cj = Yliez xihj)- 
We refer to both row and column sums as the line sums of F. We will usually only consider 
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finite sequences IZ = (ri, r2, . . . , r^) and C = (ci, C2, . . . , c„) of row and column sums tliat 
contain all the nonzero line sums. In this paper we will always assume that the line sums 
are monotone, that is ri > r2 > . . . > and Ci > C2 > . . . > c„. 

Given sequences of integers 71= (ri, r2, . . . , Tm) and C = (ci, C2, . . . , c„), we say that (JZ, C) 
is consistent if there exists a set F with row sums TZ and column sums C. Define hi = #{j : 
Cj > i} for i = 1,2, . . . ,m. Ryser's theorem [12] states that if ri > r2 > . . . > r^, the line 
sums (TZ, C) are consistent if and only if for each A; = 1, 2, . . . , m we have Y^^^i hi > Yli=i ^i- 
Note that by definition of hi we have Yll^i = Si=i ~ SI^i ^i- 

We will now define a uniquely determined neighbour corresponding to line sums {7Z,C). 
This is a uniquely determined set that is in some sense the closest to any set with those 
line sums. See also Section 4]. 

Definition 1. Let row sums TZ = (ri, r2, . . . , r^) and column sums C = (ci, C2, . . . , c„) be 
given, where n = ri > r2 > ■ ■ ■ > Vm and m = Ci > C2 > . . . > c„. Let hi = ^{j : Cj > i} for 
i = 1,2, . . . ,m. Then the column sums C\, C2, . . . , c„ and row sums h\,h2, . . . ,hm uniquely 
determine a set Fx, which we will call the uniquely determined neighbour corresponding 
to line sums iJZ, C) . 

Suppose line sums TZ= {r\,r2, ■ ■ ■ , rm) and C = (ci, C2, . . . , c„) are given, where ri > r2 > 
. . . > Tm and ci > C2 > . . . > c„. Let the uniquely determined neighbour corresponding to 
{TZ,C) have row sums &i > 62 > • • • > ^n- Then we define 

^ m 

a{TZ,C) = h\. 

i=l 

Note that a(7Z,C) is an integer, since 2a(JZ,C) is congruent to 

m m m m 

Y^{ri + hi) = ^ri + ^hi = 2^ri = Q mod 2. 

j=l i=l 4=1 i=l 

If we write di = hi — rj for all i, then because — Si^i have 

a = ^di = -"^di. 

di>0 di<0 

We can view the set F as a picture consisting of cells with zeroes and ones. Rather than 
(i, j) e F, we might say that {i,j) has value 1 or that there is a one at {i,j). Similarly, for 
^ F we sometimes say that has value zero or that there is a zero at {i,j). 

We define the boundary of F as the set consisting of all pairs of points {{i,j), {i',j')) such 
that 
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• i — i' and \j — j'\ — 1, or \i — i'\ — 1 and j — f, and 

One element of this set we call one piece of the boundary. We can partition the boundary 
into two subsets, one containing the pairs of points with i = i' and the other containing 
the pairs of points with j — f. The former set we call the vertical boundary and the 
latter set we call the horizontal boundary. Wc define the length of the (horizontal, vertical) 
boundary as the number of elements in the (horizontal, vertical) boundary. For a given set 
F we denote the length of the horizontal boundary by Lh{F) and the length of the vertical 
boundary by Ly{F). 



3 The construction 

In this section we will construct a set F2 satisfying given monotone row and column sums 
that are consistent. First we will describe one step of this construction. 

Let row sums TZ — (ri, r2, . . . , r^) and column sums C — (ci, C2, . . . , c^) be given, such that 
n — ri > r2 > . . . > rjn and m = ci > C2 > . . . > c„. Assume that those line sums are 
consistent. For i = 1,2, ... ,m define hi = ^{j : Cj > i} and di = hi — ri. For convenience 
we define r^+i = 6^+1 = c?m+i = 0. We have n = 61 > 62 > • • • > &m > &m+i- 

Let Fi be the uniquely determined neighbour corresponding to the line sums {TZ,C). Then 
Fi has row sums (fei, 62, ■ ■ ■ , bm) and column sums (ci, C2, . . . , c„). Moreover, in every column 
j the elements of Fi are exactly in the first Cj rows. 

If Tj = bi for all i, then Fi already satisfies the line sums {7l,C), and there is nothing to 
be done. Now assume that not for all i we have r^ — b^. Then there is at least one i with 
di > and one i with di < 0. Also, because of the consistency of the line sums the smallest 
i with di ^ satisfies di > 0. 

Let ii be minimal such that di^ > and let 12 be minimal such that di^ > and (ijj+i < 0. 
Let R'^ — {ii,ii + l,...,i2}. Similarly, let i^ be minimal such that di^ < and let 
be minimal such that di^ < and di^+i > 0. Such ^4 exists, since dm+i — 0. Let R~ — 
{is, is + 1, ... , i4}. Now di > for all i e and di < for all i e R~. 

If < we execute an A-step, while if \R'^\ > \R~\, we execute a B-step. We will 
now describe these two different steps. 

A-step. Let j be maximal such that Cj e R~^. Such a j exists, because as bi^+i < Ti^+i < 
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< bi^, there exists a column with sum i2. Define s — Cj — ii + 1; this is the number of rows 
i with ii < i < Cj. We will be moving the ones in the s cells ■ ■ ■ , (cj, j) to other cells. 

To determine to which cells those ones are moved, consider 23, 23 + 1, . . . , 23 + s — 1. Since 
^4 — ^3 + 1 = l-R^ I > l-R^I > s, we have ^3 + s — 1 < ^4, so {^3, + 1, . . . ,i3 + s — 1} C R~ . 
If Ti^+s-i > Ti^+s, then let I = {^3, ^3 + 1, . . . , is + s - 1}. 

Now suppose rj3+s_i = ri.^^s- Let ti be minimal such that ^3 < < ^3 + s — 1 and 
^ii — '^i:i+s~i- Let t2 be such that t2 > is + s and ri,^^s_i = Vp^ > Vf^+i- Since we have 
di^+i > 0, we have r^^+i < k^+i < < Vi^, hence t2 < i^- Let ^3 = ^2 + - ^3 - ^ + 1. 
As ^2 > is + -s, we have ^3 > ti + 1, and as ti < is + s — 1, we have ts < ^2- Now define 
I = {is, is + 1, . . . , ii - 1} U {is, is + 1, ■ ■ ■ , ^2}- We have |7| = {h - h) + (-^i + ig + s) = s. 

In both cases we have now defined a set / C i?^ with \I\ = s — Cj — ii + 1 and satisfying 
the following property: if i G / and i + 1 ^ /, then > rj+i. 

Now we move the ones from the rows i with ii < i < Cj to the rows i & I. This column 
will later be one of the columns of F2. We delete the column and change the line sums 
accordingly: define for i = 1,2, ... ,m the new row sums r^, which is equal to if there 
was no one in this row in column j, and equal to — 1 if there was a one in this row in 
column j. We have 

j — 1 for i < ii, 
j for ii < i < Cj, 
j — 1 for i e 7, 
j for i > Cj and i ^ I. 

Also let b'^ be the number of columns not equal to j with column sum at least i. We have 
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bi — 1 for i < Cj, 
bi for i > Cj. 



Note that the set F{, defined as Fi without column j, has row sums b[,b2,..., b'^. 

Wc now want to show that the new row sums are non-increasing and that they arc consis- 
tent, together with the column sums without column j, that is, that b[ > Y^'l^i r[ for 
k — 1,2, . . . ,m. 

Suppose for some i we have r[ < r^_,_^. Then we must have r[ — — 1 and r^_,_^ = r^+i, 
since > rj+i. So cither i = ii — lorie / and i + 1 ^ /. In the latter case we know 
Tj > Tj+i, hence r- > r^^^i- If on the other hand i = ii — 1, we have di = and dj+i > 0, so 
= bi > > Tj+i, hence r- > r^_,_]^. We conclude that it can never happen that < r^_,_i. 
So n — 1 = > r2 > . . . > r^. 
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Now we prove consistency. For i < ii we have di — 0, hence 

K-ri = {bi-l)-{n-l) = di = Q. 
For ii < i < Cj we have di > 0, hence 

K-r'i = {bi-l)-ri = di-l> 0. 
For Cj + 1 < i < ^3 — 1 we have di > 0, hence 

b'i - '''i = bi - n = di > 0. 
So for A; < ^3 — 1 we clearly have 

i^i^'i - r'i) > 0. 

i=l 

On the other hand, for k>iiwe have Y^\=i{hi — > because of the consistency of the 
original line sums, hence 

k / k \ f ^ \ ^ 

T.(^'.-r'i}-[J2^^-c^-[J2ri-c^^J2(^^-n)>0. 

i=l \i=l J \i=l / 1=1 

For i3 < i < i4 we have di < 0, so 

b'i - r\ = bi - ri = di < Hi ^ I, 

b'.-rl^bi-{ri-l)^di + l<0 Hie I. 
Hence for ^3 < A; < ^4 — 1 we have 

k ^4 14 

i=l 1=1 i=k+l 

This proves the consistency. 

B-step. Let j be minimal such that Cj + 1 G . Such a j exists, because as 623-1 > 
^"13-1 > f^is > bi:^, there exists a column with sum i^ — 1. Similarly to the A-step, we find a 
set / C i?"*" such that |/| = ^4 — Cj with the following property: if i ^ I and i + 1 G /, then 
n > n+i. 

Now we move the ones from the rows i with i G / to the rows i with Cj + 1 < i < 24. This 
column will later be one of the columns of F2. Wc delete the column and change the line 
sums accordingly. Analogously to above we prove that the new line sums are non-increasing 
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and consistent, and that the set F[ that we have left, is the uniquely determined neighbour 
corresponding to these new line sums. 

The procedure described above, which changes line sums (JZ, C) and their uniquely deter- 
mined neighbour Fi to new line sums {Tl',C') and their uniquely determined neighbour 
F[, we denote by (p. Since the new line sums satisfy all the necessary properties, we can 
apply if also to {TZ', C) and F[. We can repeat this until we arrive at a situation where the 
uniquely determined neighbour already satisfies the line sums. One by one we can then put 
the deleted columns back in the right position (first the column that was last deleted, then 
the one that was deleted before that, and so on, to make sure that the resulting set F2 has 
its columns in the right order). Every time we put back a column, the line sums change 
back to what they were before that instance of (p was applied. When all the columns are 
back in place, the line sums are therefore equal to (JZ, C) and the resulting set satisfies 
these line sums. This proves the following theorem. 

Theorem 1. Let row sums TZ = (ri,r2, . . . ,rm) and column sums C = (ci,C2, . . . ,c„) be 
given, where n = ri > r2 > . . . > and m = ci > C2 > ■ ■ ■ > Cn- Assume that the line 
sums are consistent. Let Fi be the uniquely determined neighbour corresponding to the line 
sums {7Z,C). If we start with Fi and repeatedly apply (f until this is no longer possible, and 
then put all the deleted columns back in the right position, then the result is a set F2 that 
satisfies the line sums {7Z,C). 

Now we show an example of this construction. Let m = 12, n = 11 and define line sums 
7^ = (11, 10, 8, 8, 8, 6, 6, 6, 3, 3, 3, 2), C = (12, 10, 7, 6, 6, 6, 6, 6, 6, 6, 3). 

We have 

(61,..., 612) = (11,11,11,10,10,10,3,2,2,2,1,1), 

(di, du) = (0, +1, +3, +2, +2, +4, -3, -4, -1, -1, -2, -1). 

We will now do the construction step by step, illustrated by Figure [T] The and di in 
every step are indicated in the figure. We start with the uniquely determined neighbour 
Fi, that is, the set with column sums C and row sums (61, ... , 612). 

Step 1. We have R+ = {2,3,4,5,6}, R = {7,8,9,10,11,12}. Since \R+\ < we 
execute an A-step. The rightmost column j with cj G R^ is column 11, with sum 3. We 
delete the ones in (2, 11) and (3, 11). We find / = {7, 8}, since rg > rg. So we add ones in 
(7, 11) and (8, 11). We then delete column 11. 

Step 2. We have R+ = {3,4,5,6} and R- = {7,8,9,10,11,12}. Since \R+\ < \R-\, we 
execute an A-step. The rightmost column j with Cj G i?"*" is column 10, with sum 6. 
We delete the ones in this column in rows 3, 4, 5 and 6. Since rio = rn, we cannot use 
/ = {7, 8, 9, 10}. Instead we take / = {7, 8, 10, 11}. This works since rg > rg and rn > r^- 
So we add ones in column 10 in rows 7, 8, 10 and 11. We then delete column 10. 
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(a) Step 1. 
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(b) Step 2. 
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(c) Step 3. 



(d) Step 4. 
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(e) Step 5. 
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(f) The set F2 with its boundary. 



Figure 1: The construction of the set F2. The ones are indicated by white circles. The dashed 
circles are ones that are deleted in that step, while the black circles are ones that are newly added 
in that step. The numbers directly next to each ^gure are the row sums, while the numbers next 
to that are the dj. 



Step 3. In row 10, the new row sum is 2, while the new bio is also 2. So the new diQ is 0. 
This means that while is still equal to {3,4,5,6}, we now have = {7,8,9}. Hence 
l-R"*"! > \R~\ and therefore we execute a B-step. The leftmost column j with Cj + 1 G R~ 
is column 3 with sum 7. So we add ones in (8,3) and (9,3). As = = r^, we cannot 
take I = {6,5}, but we have to take / = {6,3}. Hence we delete ones in (3,3) and (6,3). 
We then delete column 3. 

Step 4. We have i?+ = {4,5,6} and R~ = {7,8}. As > \R'\, we execute a B-step. 
The leftmost column j with Cj + 1 E R~ is column 3 (which was originally column 4) with 
sum 6. We add ones in (7, 3) and (8, 3). As = r4^, we take / = {6, 4}, so we delete ones 
from (4,3) and (6,3). We then delete column 3. 

Step 5. We have R^ = {5,6} and R^ = {11, 12}. As < \R^\, we execute an A-step. 
The rightmost column j with cj G -R"*" is column 7 (which was originally column 9) with 
sum 6. We deletes ones from (5, 7) and (6, 7), and we add ones in (11, 7) and (12, 7). We 
then delete column 7. 

Now all di have become 0, so we are done. We put back the deleted columns in their original 



places and find the set F2 that satisfies the original line sums, see Figure 1(f) 



4 Boundary length of the constructed solution 

In this section we prove upper bounds on the length of the boundary of the set that results 
from the construction described in the previous section. 

Theorem 2. Let row sums TZ = (ri,r2, . . . ,rm) and column sums C = (ci,C2, . . . ,c„) be 
given, where n = ri > r2 > . . . > and m = ci > C2 > ■ ■ ■ > Cn- Assume that the line 
sums are consistent. Let a = a(JZ,C). For the set F2 constructed in Theorem^ we have 

Lh{F2) <2n + 2a, ^^(Fa) < 2m + 2a. 



Proof. Let Fi be the uniquely determined neighbour corresponding to the line sums {TZ, C). 
Starting with Fi, we apply (p repeatedly, moving ones in several columns accordingly and 
deleting those columns. After that, to analyse what happens to the boundary, we start 
again with Fi and repeat the entire procedure, moving exactly the same ones, but this 
time keeping the columns that were supposed to be deleted. 

The length of the horizontal boundary of Fi is equal to 2n, since there are n columns that 
each contain one connected set of ones. The length of the vertical boundary of Fi is 2m. 
Note that the ones that are moved when applying (p are always deleted from a row i with 
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di > and added to a row i with di < 0. In fact for each row i with di > ones are 
deleted exactly di times during the construction, and for each row i with di < ones are 
added exactly —di times. Therefore the total number of ones that are moved is equal to a. 
We now want to show that when in one application of (p exactly s ones are moved, both 
the horizontal and vertical boundary do not increase with more than 2s. Prom this the 
theorem follows. 

We will only consider what happens at an A-step; the other case is analogous. So suppose 
we execute an A-step and move s ones, while either the horizontal or vertical boundary 
increases by more than 2s. First consider the horizontal boundary. Since the ones in the 
rows i with ii < i < cj are removed, and there never was a one in {cj + this does 
not yield any additional boundary. Adding the ones in the rows i with i E I may yield 
additional boundary, but only 2 for each one that is added, so at most 2s in total. 

So we may assume that the vertical boundary has increased by more than 2s. Adding the 
ones leads to additional vertical boundary of at most 2s, so deleting the ones must also 
have led to additional boundary. This means that there was a one in which is now 

deleted, while there are still ones in {i.j — 1) and + 1). As di > 0, those ones cannot 
have been added during earlier steps in the construction, so they must have been there 
from the beginning. This means in particular that Cj+i > i > ii, while also Cj+i < cj < i2, 
so Cj+i e R'^. But j was chosen maximally such that cj e R'^, so apparently column j + 1 
was in the original construction deleted in an earlier application of (p. 

Suppose this earlier application has been an A-step. Since rows / with < at some 
point in the construction can never have > at a later point in the construction, we 
know that all rows I with ii < I < Cj+i were contained in in this earlier application of 
(f. In particular should the one in + 1) have been moved during this step. So this is 
impossible. 

Now suppose that the earlier application has been a B-step. Then column j + 1 can only 
have been chosen to execute this step in if (icj_,_^+i < 0. Since Cj+i < cj and dc^ > (now, 
and therefore also earlier), we then must have Cj = Cj+i. Hence rfcj+i < 0, which means 
that to execute this B-step column j, rather than column j + 1, should have been chosen. 
So this case is impossible as well. 

We conclude that the vertical boundary has increased by at most 2s as well, and this 
completes the proof of the theorem. □ 

In light of this theorem it is interesting to note that a cannot become arbitrary large while 
n and m are fixed. In fact, we have the following result. 

Proposition 1. Let row sums TZ — (ri, r2, . . . , r^) and column sums C = (ci, C2, . . . , c„) 
be given, where n — ri > r2 > ■ ■ ■ > rm and m = ci > C2 > . . . > c„. Assume that the line 
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sums are consistent. Let a = a{TZ,C). Then 

(m — l){n — 1) 



a < 



4 



Proof. For i = 1,2, ... ,m, let bi = #{j : Cj > i} and di = bi — ri. Let a be the number of 
rows (indices i) with di > and b the number of rows with di < 0. We assume a > 0, so 
a,b > 0. Define d~^ = max{di : di > 0} and d~ = max{— c?j : di < 0}. We have bi = n = ri, 
so di = 0, hence a + b < m — 1. 

Now we prove that d'^ + d~ < n — 1. Let k and / be such that bk — r^ = d^ and ri — bi = d~ . 
First suppose k < I. Then since ri > r2 > . . . > and &i > &2 ^ • • • > we have 
bi>bk = bk - Tk + Vk = d^ + Vk and > -6/ = ri - bi - ri = d' - ri, hence 

(i+ + < (6i - Tk) + (-6m + Ti) <bi-hm<n-l. 

If on the other hand k > I, then ri > ri = ri — bi + bi = d^ + bi and — > — Tfc = 
bk — f^k — bk = d'^ — bk, and hence 

d^ + d' < {-Vm + bk) + (ri -bi) <ri-rm<n-l. 



Now note that we have 

a = ^d, = ^{-di). 



di>0 di<0 



SO 



a' 



\ ) ^ • ■ = {a-b) {d^ ■ d-) 

a + b\'^ f d^ + d^ \ /m — lA^/n — 1^^ 



< — < 



Therefore 



a < 



2 J - \ 2 
(m — l)(n — 1) 



□ 



In case of large a, the construction of Theorem [T] actually gives a much smaller horizontal 
boundary than the bound in Theorem [2| as the following theorem shows. 

Theorem 3. Let row sums TZ = {ri,r2, . . . ,rm) and column sums C = (ci,C2, . . . ,c„) be 
given, where n > 2, n = ri > r2 > . . . > Vm and m = ci > C2 > . . . > c^- Assume that the 
line sums are consistent. For the set F2 constructed in Theorem\l\we have 

Lh{F2) < 4n - 4. 



11 



Proof. We will prove this by induction on n. Let a — a{TZ,C). If a > 0, then there are h 
and I2 such that 2 < li < I2 and di^ > and di^ < 0. Then 

bi>bi,>ri, + l>ri, + l>bi, + 2> 1 + 2^3. 

Hence n > 3. So when n = 2, we have a — and the construction gives F2 — Fi, with 
Lh ^ 2n = An - 2n = An - 4. 

Now let A; > 3 and suppose that we have proved the theorem in case n < k. Let n — k. Let 
Fi be the uniquely determined neighbour corresponding to the line sums {TZ, C) . We apply 
(f to Fi once. Assume without loss of generality that an A-step is executed in column j. 

First suppose that / consists of consecutive numbers. Then after moving the ones in column 
j, the length of the horizontal boundary in this column is equal to 4. When we delete this 
column, we are left with A; — 1 columns, so we can apply the induction hypothesis, which 
yields that the total length of the horizontal boundary at the end of the construction will 
be 

< 4(A; - 1) - 4 + 4 = 4A; - 4. 

Now suppose that / does not consist of consecutive numbers. Then we know that / is of 
the form I — {is, is + 1, . . . , — 1} U {ts, is + 1, ... , ^2}- So after moving the ones, the 
length of the horizontal boundary in column j is equal to 6. Also, we know in particular 
that the one in {cj,j) was deleted and a one was added in (is, j)- 

The new parameters, after moving the ones and deleting column j, we denote by r^, b'^ and 
d'^. The construction will in later steps execute an A-step in at most d'^^_i columns with 
sum is — 1 and a B-step in at most —d[^ columns with sum is — 1. On the other hand, we 
currently have b\^_i — b\^ columns with sum is — 1. 

We know that Vi^-i > ri^, and r^^ = rjg — 1. Both in the case cj = is — 1 and in the case 
Cj < is — 1, we have r-^.^ = ri^-i, so 

(6:3_i - 6U - (4_i - 4) = r:3_i - r; = n3_i - r,3 + 1 > 1. 

This means that there is at least one column with sum is — 1 in which none of the later 
steps of the construction will be executed. This column will at the end of the construction 
therefore still have a horizontal boundary of length 2. If we delete this column entirely 
and then do the construction, exactly the same steps will be carried out. After all, the 
deleted column would never have been chosen to execute a step in anyway; also, deleting 
the column does not influence the choice of the set / in each step of the construction, as the 
only difference between the row sums of two consecutive rows that is changed, is between 
rows is — 1 and is, but as di^^i > and di^ < 0, these rows will never both be in or 
both be in R~ . 
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By applying the induction hypothesis to the new situation with n = k — 2, we find that 
the total length of the horizontal boundary at the end of the construction will be 



Lft < 4(A; - 2) - 4 + 6 + 2 = 4A; - 4. 
This completes the induction step. □ 

Unfortunately, we cannot prove a similar result for the vertical boundary. In fact, we can 
find examples for which our construction gives a vertical boundary with a length as large 
as + |m + ^, see Example [S However, we believe that there always exists a solution 
with a small boundary length, both horizontal and vertical. 

Conjecture 1. Let row sums TZ = (ri, r2, . . . , r^) and column sums C = (ci, C2, . . . , c„) be 
given, where n = ri > r2 > . . . > and m = ci > C2 > . . . > c„. Assume that the line 
sums are consistent. There exists a set F3 with line sums {TZ, C) for which 

LhiF3)<4n-4, L,(F3)<4m-4. 



5 Examples 

We give two families of examples for which we can prove that the construction of Theorem 
[1] gives the smallest possible length of the boundary. 

Example 1. Let the number of columns n be odd and let m = n. Define line sums 
C = 7^ = (n, n - 1, n - 1, n - 3, n - 3, . . . , 4, 4, 2, 2). 

We calculate 

(61, 62, • • • , &n) = (n, n, n - 2, n - 2, . . . , 3, 3, 1), 

(rfi, c?2, . . . , dn) = (0, +1, -1, +1, -1, . . . , +1, -1, +1, -1). 

So a = a(Jl,C) = Theorem [2] tells us that the set F2 constructed with Theorem [l] 
satisfies 

Lh{F2) <2n + 2a = 3n-l, ^^(Fs) < 2m + 2a = 3n - 1. 

On the other hand, by ^ Corollary 1] we know that for any set F with these line sums, 
we have 

Ti — 1 

Lh{F) >2n+ ^— ■ (1 - (-1)) + 2 ■ = 3n - 1, 

and by symmetry also Ly{F) > 3n — 1. This shows that F2 has the smallest boundary 
among all sets F with these line sums. See for the constructed set F2 in the case that n = 9 



Figure 2(a) (This example is in fact a slightly modified version of [6l Example 3].) 
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(a) Example [l] with n = 9. 
The horizontal and vertical 
boundary both have length 
26. 
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(b) Example [2] with n = 9. 
The horizontal and vertical 
boundary both have length 
32. 
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(c) Example [5] with fc = 3. 
The vertical boundary has 
length 50. 



Figure 2: The constructed sets F2 for some of the examples. 



Example 2. Let m = n > 2. Define line sums 

C = 7^ = (r^,2,2,2,...,2,2). 

We calculate 

(61,62,- •• ,&n) = 
(di, d2,..., dn) = (0, n - 2, -1, -1, . . . , -1). 

So a = a{TZ,C) = n — 2. Theorem [2] tells us that the set F2 constructed with Theorem [l] 
satisfies 

Lh{F2) < 2n + 2a = 4n - 4, ^^(Fs) < 2m + 2a = 4n - 4. 

On the other hand, by [6, Corollary 1] we know that for any set F with these line sums, 
we have 

Lh{F) >2n + 2{n - 2) = 4n - 4, 

and by symmetry also L^{F) > An — A. This shows that F2 has the smallest boundary 
among all sets F with these line sums. See for the constructed set F2 in the case that n = 9 
Figure 2(b) (This example is in fact a slightly modified version of Example 4].) 



We can generalise Example |2] to larger a, in which case the bound of Theorem [2] is no 
longer sharp. However, in this case we can use Theorem [3] to prove that the horizontal 
boundary is the smallest possible, as shown below. 

Example 3. Let k he a. positive integer and let m = kn — k + 1. Define line sums 
C = {kn- k + l,k + l,k + l,...,k + l,k + l), 7^ = (n, 2, 2, . . . , 2). 
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We calculate 

{bi,b2,...,bm) = ( n,n,^..,n , l,l,^..,l ), 

fc+l times kn—2k times 

(c?i,c?2, . . . = (0,n - 2,n - 2, . . . ,n - 2, -1, -1, . . . , -1). 

V ^ V 

fc times kn—2k times 

Theorem |3] tells us that the set F2 constructed with Theorem [T] satisfies 

Lh{F2) < 4n - 4. 

On the other hand, by [6, Corollary 1] we know that for any set F with these line sums, 
we have 

LhiF) >2n + 2{n - 2) = 4n - 4. 

This shows that F2 has the smallest horizontal boundary among all sets F with these line 
sums. 

The next example shows that the upper bound on a given in Proposition [T] can be achieved. 
Example 4. Let A; be a positive integer and let m = n = 2/c + 1. Define line sums 

C = n= {2k + l,k + l,k + l,...,k + l). 

We calculate 

(6i,62,...,U = ( 2A: + l,2fc + l,...,2fc + l ,l,l,...,l), 

V '' V ' 

fc+1 times fc times 

{di,d2, ...,dm) = {0, k, k, . ^.., k , -k, -fc^. . . , -k) . 

k times k times 

Hence 

,9 (m — l)(n — 1) 
a = a{n,C) = P = ^ 2 ■ 

Finally we show by an example that the vertical boundary of the set F2 constructed in 
Theorem [1] can become quite large, so it is not possible to prove a similar result as Theorem 
|3]for the vertical boundary. 

Example 5. Let A; be a positive integer and let m = 3k + 1, n = 3k. Define line sums 
C = (3k + l,k + l,k + l,...,k + l), 7^ = (3k, k + l,k + 1, . . . ,k + l,k,k, . . . ,k). 

^ V ' ^ v ' 

2k times fc times 

We calculate 

(61, &2, • • • , bm) = ( 3/c, 3k , 3k , 1,1,..., 1), 

""^ V ' ^ v ' 

fc+1 times 2fc times 
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{di, d2, . . . , dm) — 

(0, 2k-l,2k-l,...,2k-l, -k, -k,..., -k, -{k - 1), -{k - 1),. . . ,-{k - 1)). 

k times k times k times 

Hence a = a{n, C) = 2k^ - k. 

The construction executes 2k — 1 times an A-step, in each of the columns 3k, 3k — 1, . . . , 
k + 2. In the first step (and every odd-numbered step after that) we have I = {k + 2, k + 
3, ... ,2k + 1}. At the beginning of the second step, however, the row sums in rows k + 2, 
k + 3, . . . , 3k + l are all equal, so we have / = {2k + 2, 2A; + 3, . . . , 3A; + 1}. The same holds 
for every other even-numbered step. This means that at the end of the construction, the 
vertical boundary in each of the rows k + 2, k + 3, . . . , 2k + 1 will be equal to 2(fc + 1), 
while the vertical boundary in each of the rows 2k + 2, 2k + 3, . . . , 3k + 1 will be equal to 
2k. Adding the boundary of 2 in each of the rows 1, 2, . . . , A; + 1, we find 

L^{F2) = {k + 1) ■ 2 + k ■ 2{k + 1) + k ■ 2k = Ak^ + ik + 2. 

This is not linear in m = 3/c + l. It is in fact equal to |m^-|- |m+ ^. See for the constructed 
set F2 in the case that k = 3 Figure |2(c) 

It is clear that in fact there exists a set F with the same line sums, but with a much smaller 
vertical boundary, which supports Conjecture [1} 



6 Generalising the results for arbitrary ci and ri 

In all results of the previous sections, we used the condition that ci = m and ri = n. 
This is purely for convenience; it is not a necessary condition. We can easily generalise the 
results for the case that these conditions do not necessarily hold. 

Consider a set F with row sums 7^ = (ri, r2, . . . , r^) and column sums C = (ci, C2, . . . , c„), 
where ri > r2 > . . . > and Ci > C2 > . . . > c„, but not necessarily ci = m and ri = n. 
Let F' be a set that is equal to F, except that we add a full row with index and a full 
column with index 0, i.e. 

F' = FU {(0,j) : < J < n} U {(i,0) : 1 < i < m}. 

The row sums of F' are 

= {r'o,r[,r'2,. . .,r'^) = (n, ri + 1, r2 + 1, . . . , + 1). 

and the column sums of F' are 

C = (c'o, C[, C2, . . . , C'„) = (m, Ci + l,C2 + l,...,Cn + l). 
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It is easy to see that a{JZ',C') = a(7l,C). Now consider the length of the horizontal 
boundary. For every j with (1, j) G F, the horizontal boundary in column j of F' is equal 
to the horizontal boundary of column j in F. For every j with (1, j) ^ F, however, the 
horizontal boundary in column j of F' is 2 larger than the horizontal boundary in column 
j of F. (This also holds for column 0, where the horizontal boundary of F had length 
and the horizontal boundary of F' has length 2.) Hence 

L,(F')=L,(F) + 2(n + l-ri). 

Analogously, we have 

L,{F')=L,{F) = 2{m + l-c,). 

By applying Theorems [2] and [s] as well as Proposition [l] to F' (with n + 1 columns and 
m + 1 rows), we acquire the following results. 

Proposition 2. Let row sums TZ = (ri, r2, . . . , r^) and column sums C = (ci, C2, . . . , c„) 
be given, where ri > r2 > ■ ■ ■ > and Ci > C2 > . . . > c„. Assume that the line sums are 
consistent. Let a = a{TZ,C). Then 

mn 

Theorem 4. Let row sums TZ = (ri, r2, . . . , r^) and column sums C = (ci, C2, . . . , c„) be 
given, where ri > r2 > . . . > and ci > C2 > . . . > c„. Assume that the line sums are 
consistent. Let a = a{TZ,C). Then there exists a set F2 satisfying these line sums such that 

Lh{F2) < min( 2ri + 2a, 2ri + 2n - 2 ) 

and 

K{F2) <2ci + 2a. 



References 

[1] E. Balogh, A. Kuba, C. Devenyi, A. Del Lungo, Comparison of algorithms for re- 
constructing hv-convex discrete sets. Linear Algebra and its Applications 339 (2001) 
23-35. 

[2] E. Barcucci, A. Del Lungo, M. Nivat, R. Pinzani, Reconstructing convex polyominoes 
from horizontal and vertical projections. Theoretical Computer Science 155 (1996) 
321-347. 

[3] M. Chrobak, C. Diirr, Reconstructing hv-convex polyominoes from orthogonal projec- 
tions. Information Processing Letters 69 (1999) 283-289. 

[4] G. Dahl, T. Flatberg, Optimization and reconstruction of hv-convex (0,l)-matrices, 
Discrete Applied Mathematics 151 (2005) 93-105. 



17 



[5] B.E. van Dalen, On the difference between solutions of discrete tomography problems, 
Journal of Combinatorics and Number Theory 1 (2009) 15-29. 

[6] B.E. van Dalen, Boundary length of reconstructions in discrete tomography, 
|arXiv:1006.4449 | [math. CO] (2010) 25 pp. 

[7] R.J. Gardner, P. Gritzmann, D. Prangenberg, On the computational complexity of 
reconstructing lattice sets from their X-rays, Discrete Mathematics 202 (1999) 45-71. 

[8] S.B. Gray, Local properties of binary images in two dimensions, IEEE Transactions 
on Computers 20 (1971) 551-561. 

[9] G.T. Herman, A. Kuba, editors. Discrete Tomography: Foundations, Algorithms and 
Applications, Birkhauser, Boston (1999). 

[10] G.T. Herman, A. Kuba, editors. Advances in Discrete Tomography and Its Applica- 
tions, Birkhauser, Boston (2007). 

[11] A. Rosenfeld, Connectivity in digital pictures. Journal of the Association for Comput- 
ing Machinery 17 (1970) 146-160. 

[12] H.J. Ryser, Combinatorial properties of matrices of zeros and ones, Canadian Journal 
of Mathematics 9 (1957) 371-377. 

[13] G.J. Woeginger, The reconstruction of polyominoes from their orthogonal projections. 
Information Processing Letters 77 (2001) 225-229. 



18 



